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ABSTRACT 

The  equations  of  motion  are  developed  and  solved  numerically 

for  the  trajectory  of  a  spherical  particle  passing  through  a  Prandtl- 

Meyer  expansion  fan.     The  effect  of  a  change  in    v    and    n    is  shown, 

where    v    is  the  ratio  of  specific  heats  of  the  gas,   and    n/2    is  the 

exponent  in 

n/2 

o  ~o 

the  as  sunned  viscosity  temperature  relationship. 

It  is  demonstrated  that  for  particles  the  order  of  a  micron  in 
diameter,  slip  flow  will  exist,  and  a  proposed  correction  to  the  drag 
to  account  for  this  discontinuous  nature  of  the  flow  is  investigated. 

The  results  are  plotted  showing  particle  trajectory  profiles 
and  the  components  of  the  relative  velocity  or  slippage  velocity  of 
the  particle  for  flow  deflection  angles  up  to  70  degrees. 
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I.    INTRODUCTION 

Heterogeneous  flow  of  a  solid  particle  -  gas  mixture  has  been 
a  topic  of  increased  interest  in  recent  years.     The  addition  of  light 
metals  to  solid  rocket  fuels  has  focused  attention  on  the  non-equilibri- 
um nature  o£  two-phase  flows  in  expansion  processes,   and  numerous 
papers  have  been  written  regarding  this  phenomenon. 

Early  studies  of  two-phase  flov/3  were  involved  with  investi- 
gations of  water  droplet  trajectories  around  a  variety  of  body  shapes 
in  connection  with  icing  of  aircraft  surfaces.     More  recently,   hetero- 
geneous  nozzle  flow  iias  been  examined.     Gilbert,    Davi3,   and  Altman 
(reference  1)  related  rocket  thrust  losses  to  solid  particle  sizes  and 
solved  linearized  nozzle  equations  assuming  Stokes  drag  law.     Kliegel 
(reference  2)  treated  simultaneous  velocity  and  thermal  lag  along  with 
the  effect  of  particle  lag  on  gas  properties  in  a  one -dimensional  noz- 
zle.    Rannie  (reference  3)  showed  that  for  most  cases  of  interest, 
results  of  equal  accuracy  to  numerical  integration  can  be  obtained 
through  a  perturbation  analysis  of  temperature  and  velocity  lags  in  a 
one-dimensional  nozzle.     Rannie  also  introduced  slip  flow  corrections 
for  drag  and  heat  transfer.     An  excellent  summary  of  recent  advances 
in  heterogeneous  nozzle  flow  was  presented  at  the  American  Rocket 
Society,    Solid  Irropellant  Rocket  Conference  in  January  1962,  by  R. 
F.   Hoglund  (reference  4). 

Several  problems  relating  to  supersonic  flow  of  two-phase 
systems  have  also  been  examined.     Carrier  (reference  5)  investigated 
a  solution  to  particle  flow  through  normal  shocks,   and  oblique  shocks 
were  treated  by  Morgenthaler  (reference  6). 
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The  behavior  of  a  single  particle  passing  through  a  Prandtl- 
Meyer  expansion  furnishes  an  example  of  heterogeneous  flow  that  can 
be  analysed  in  a  straightforward  manner,   and  serves  to  illustrate 
some  of  the  problems  and  important  parameters  of  two-phase  flow  in 
general.     This  paper  presents  the  development  of  the  equations  of  mo- 
tion of  small  spherical  particles  flowing  with  a  gas  that  experiences  a 
Prandtl-Meyer  expansion,   and  the  numerical  solution  of  these  equa- 
tions.    The  effects  of  the  temperature  dependence  of  viscosity,   changes 
in  the  ratio  of  gas  specific  heats,   y  ,   and  slip  flow  correction  to  the 
drag  coefficient  were  investigated  by  varying  these  factors  in  succes- 
sive solutions.     The  results  of  this  investigation  have  direct  bearing 
on  the  exhaust  of  an  underexpanded  rocket  nozzle  as  the  flow  expands 
around  the  lip  of  the  nozzle. 


II.    DEVELOPMENT  OF  BASIC  EQUATIONS 

An  ideal  two-dimensional  Prandtl-Meyer  expansion  about  a 
sharp  corner  in  a  semi-infinite  channel  is  the  assumed  model  for 
this  investigation.     A  cylindrical  coordinate  system  with  origin  at 
the  vertex  of  the  expansion  fan  is  used  to  describe  the  position  of  a 
spherical  particle  passing  through  the  fan.     Denoting  the  radial  dis- 
tance from  the  origin  to  the  center  of  the  particle  as  ,    and  the 
angular  displacement  of  the  particle  from  the  initial  characteristic 
a3          (Pig-    1)»   the  conservation  of  momentum  equations  can  be 
written 

(1) 


_L(rr?r20)  -    -f   GTT  cr^u  (t^  -  re  )  r        .  (2) 

It  is  assumed  that  the  only  force  acting  on  the  particle  is  a  drag  force 
which  can  be  represented  by  Stokes  drag  law  modified  by  the  factor 
^    .      The  radial  and  tangential  components  of  the  gas  velocity  are      u*, 
and     V^     respectively,   and  M    is  the  local  value  of  the  gas  viscosity. 
The  particle  has  mass    m    and  radius    cr  . 

It  is  convenient  to  rewrite  the  equations  with    ©■  the  independ- 
ent variable.     Choosing       i_  -    m  r2-Q-  i     the  angular  momentum,   and 
Yr      as  the  desired  dependent  variables,   the  time  derivative  can  be 
written 

dt  mrz  d-8-         J    • 

Equation  1  then  becomes 


-4- 

2 


f  (orrery  (^yi) 


6r 
6  9 


Putting      —   =  -  r*  =4        ■         a«d  multiplying  by  m  r* 

da  d  0 

I  2   ci^'/r)      .  d('/r)  (  d_L  \     +     LZ 

L      del  ds     lde'        r 

+-  {(oVfcr/umr^j  -hfGTrtr/xLr2  c!l!±l  =  0 


(3) 


Equation  2  is  simply 

ik   =f67rcru^-^    _  r^  )  <4> 

d  &  '    v       L  / 

Replacing    d  L       in  the  second  term  of  equation  3  by  equation  4  we 
have 


4-f  6TTCTA77  r2C^5    =   O  (5) 

Non-dimensional  parameters  are  defined  as 

R    —      ~J~  and  A  =     r  (6) 

)  rn  c  ^  A 

/y  _  i  " 
,    y     is  the  ratio  of  gas  specific  heats,   and    C 

is  the  theoretical  maximum  expansion  velocity  determined  by  the 

equation 

C2  =    2cpT0  (if  ^y~m02)    .  (7) 

The  characteristic  length,  'f    ,   has  the  value 


*5  = 


m  c  X2  (8) 


GTTCTfJ, 
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The  subscript  zero  indicates  initial  conditions  of  the  gas  or  particle 
as  they  enter  the  expansion  fan. 

Introducing  equations  6  into  equation  5,   we  obtain 


i[^k}^f^^iR^iw-^]-o  , 


V    do 


and  dividing  by       ^>  z  &  Tf  (T  m  jU  < 


c    , 


L    d£2        Rj      r      P*L      c      d#     +    ~£- 


=.  o 


(9) 


Equation  4  can  be  written  in  terms  of  the  dimensionless  parameters 


as 


m 


R         'XT, 


cyA^  -ffitroyi^Ofr^ 


-0 


or 


(10) 


The  conservation  equations  of  frictionless  flow  yield       jfj    , 

c 

_    3.      ,   and     -l—       as  functions  of    £■  .      The  equations  may  be  writ- 

C  71 

ten  in  terms  of  our  assumed  coordinate  system  as 


Continuity:  p  (u  3  4-  1^2  J  f     -^    ^    = 


Momentum:      H 


de  J         J  de 

2 


0 


^j.  cl  u.a    _     3ZL   - 


r    de- 


o 


r 


ore 


r  V     D      d  e 


of   d£ 

fr  da 


(ii) 

(12) 
(13) 
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Energy:        a2  =  * /?3  T    -    *=L  ( Cz  -  U*  -  V)     .  (14) 

Equations  II  and  13  can  be  recombined  to  eliminate  the  term  in  the 
parentheses  and  show  that 

1/g    =    CL     .  (15) 

Using  this  result  in  equation  14,       Vj    can  be  found  as  a  function  of 
^-3       which  can,   in  turn,   be  used  in  equation  12  to  give 

^    =:    sin(AO  4-o()        < 

(16) 

^-  -  Xcos(xe  *-*)     . 
c  v  y 

The  value  of  the  constant    &(    can  be  found  by  evaluating  equations  16 
at  the  initial  conditions    Q  —  O      t    1/^  =   cxa       t   antj  using  equation 
14  to  give 

COS^=     X>(U*->!    -r     I  '  (l?) 


The  quantity      r=~       can  be  found  by  applying  equation  14  to  show  that 


71    '-r<tl-m 


or  by  equations  16 

T   ___     cosz(xe  -t-oC) 
T°  Cos2  (X 


(13) 


It  is  assumed  that  viscosity  varies  with  temperature  raised  to  a  frac- 
tional power,   hence 


£=  [J0)z 


{19) 


-7- 
or  from  equation  18 

A  COS  "l     XO     roi) 

,uo  COS"  oi 

where     n    will  normally  be  between  1.  0  and  2.  0. 

The  equationsof  particle  motion  (9  and  10)  can  now  be  written 
entirely  in  terms  of  the  dimensionless  variables    R   ,     ^    ,     &    ,   the 

drag  correction  factor    -f    ,   and  the  constants   < -'■<   and      A   .     Letting 

i  /      \ 
(      y  denote      2-i they  become 

d  e 

+    SIP  (,\0  t-o(  )      -:     o 


(21) 


At   • 


acosn(Ae-Hoc) 

cos*  <x 


_-"  A   COS(Ae+o()  -  Aj  =.  O    .    (22) 


A  value  of  £    for  spheres  in  the  subsonic  slip  flow  regime  has 
been  proposed  by  \V.    D.   Rannie  (ref.  3  )  which,   for  relative  Mach  num- 
bers less  than  unity  and  for  Reynolds  numbers  less  than  about  0.  5, 
takes  the  form 

L   ^     ,.o  +   4-OSy^     —      i  (23) 

where     Ms      and     /?e     are  the  relative  >'ach  number  and  Reynolds 
num.be r  respectively. 

,   which  is  related  to  the  ratio  of  molecular  mean  free 

path  to  the  sphere  diameter, can  be  evaluated  as  follows: 
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J^l  =        >M*        UL)(-2i)(& 
Re  '"  Zaaf0  r    (yu0A  a/lp 

From  equations  15  and  16 

Q6  COS^ 


{24} 


a.  cosO  e  +  <x) 

Isentropic  relationship  between  density  and  temperature  along  with 
equation  19  gives 


(25) 


2 


ft 


COS  <x 


Cos(A  ©  +  <*) 


Y 


Using  equations  20,    25,   and  26  in  equation  24,   we  fjet 


R 


//<= 


cos°( 


,-^-n 


,?  a,  a8o-[cos(XOt  q()_ 


(26) 


(27) 


and  equation  23  becomes 


—  —   I  + 
■f 


zozJT^\     cos* 


ou  <r 


cos(X©  +  o() 


f>-n 


(28) 


The  drag  correction  factor  cannot  be  written  in  terms  of  the  dimen- 
sionless  similarity  parameters    R    and    A     ,   but  is  dependent  on  the 
particle  size  and  initial  conditions. 

It  is  convenient  to  relate  the  trajectory  of  the  particle  to  that 
of  a  particle  having  zero  relative  velocity,  i.  e. ,  moving  alone;  a  gas 
streamline.     Designating  this  path  as       r"q    =     r^  (  0  )    }    we  can  write 


d-t  ^        do 


de 
dt 


and  since 


de  _     % 


d-b 

de 
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v/e  have 


r3 


or. 


3 


Integrating  this  equation  we  obtain 


(29) 


r. 


Ro 


COS      <* 


From  the  definitions  of     f^    and     A    it  follows  that 


3 


A 


vhich,   when  combined  with  equations  16  and  30,    gives 


(30) 


A 


COS    <* 


cos(xe  +  o() 


(31) 


Two  new  dimensionless  parameters  are  defined 

x  =      R' 


R 


y  =    A 

A, 


(32) 


from  which  it  follows  that 


R 


_X 


"  cosfA 

e  +  o()~ 

COS 

°C 

xz 


(33) 


(r)   ~  Re 


,   [cos(xe-Kx)"|M      / 


COSo<. 


x  -  ^  tan(A©+o() 

A 


»      / 


(34) 
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R 


R»        Cos  <* 


»/v- 


X 


tan  (\0+e() 


+  x{(^-Ot^W00-,]|      ■ 


(35) 


A     -    YA 


COS    oi 


c  o s (x  e  +  ©c ). 


A2 


(36) 


/  [     COSo( 

A   =   A0 


.cos(xe  +  o() 


x* 


V  +  y(-r  -X)tan(Xe+cx) 


L. 


(37) 


Also  from  equations  6  and  16 

A0    -    R0 cos    ^ 


pi ) 


Writing  equations  21  and  22  in  terms  of  X   ,  /   parameters 
using  equations  33  through  38,   we  obtain  the  result 

(cos<*)*£~n 


^ARe(^ 


yv(co5(AeM))>-h  + 


(39) 
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y'+     V  (4-~  X)  tan(A©  +  c<) 


-fARA  x3Y   /(cosCXe  +  oi))^ 


(40) 


It  is  assumed  that  the  particle  is  moving  with  no  relative  velocity  prior 
to  reaching  the  first  characteristic  line  of  the  expansion.     The  initial 

conditions  are  then 

X    (0)      zr         I.O      , 

y  (o)  =     i  o   ,  (41) 

tf'(O)    =        o      • 

Equations  28,    39,   40,    and  41  constitute  the  complete  statement  of 
the  problem. 
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m.    AUXILIARY  EQUATIONS 

In  order  to  specify  the  trajectory,  the  velocity  as  well  as  the 
position  of  the  particle  should  be  known.  The  important  parameters 
in  terms  of    X   and    y  are  then 


1L  -    A  .   R.%    -  x  y 


-UJ 


R         A 


(42) 


3.  -  £L  .& 


i  r  cos  <*    -i 


*  L  CoS(*6+oO. 


(43) 


which,   by  33,    34,   and  16,   is 

-  =  xyfi  --7  x  cot(Ae  +  o() 


(44) 


The  relative  velocity  or  3lip  velocity  of  tiie  particle  relative 
to  the  absolute  gas  velocity  is  simply 


a3  -  u 


(45) 


The  relative  Mach  number  is  then 


Ms  = 


a 


-t- 


Us 

0. 
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or 


/Mc   — 


2H 
i/*» 


U«  \2 


VlTs 


(46) 
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IV.   .SOLUTION  OF  EOUATIONS 

Since  equations  39  and  40  are  non-linear,   a  numerical  solu- 
tion is  indicated  for  arbitrary  values  of     R0    .     For  either    R0  >>  l.o 
or    R0<<    '°      i   the  possibility  of  a  small  perturbation  solution 
exists;  however,   this  possibility  was  not  investigated. 

Equations  21  and  22  were  first  solved  utilizing  a  desk  calcu- 
lator and  a  simple  iteration  procedure  where 


(i)A 


R) 


W  -L 


■+■ 


// 


a, 


&)',=  (*" 


+ 


II 


(47) 


A2  -^    At  +    ^1  A       + 


"1 


r2        *    ^ 

\S  -   A,    4- A,  Aej_ 


A6 


For  these  calculations  it  was  assumed  that    £  - ).  o     (Stokes  drag  law 
applies),     1$  -1-4    ,      n  =  20      ,   and     Mo   -   y/27       .      Various  values 
of      R0       between  0.  1  and  100  were  investigated;  however,   it  was 
found  that  the  iteration  interval     A  B  ;  had  to  be  decreased  as      fR 
increased  to  keep  the  process  convergent.     For  Ra  very  small,   the 
iteration  interval  had  to  be  made  very  small  again  to  get  meaningful 
result  (proper  divergence  from  the  zero  deflection  line). 
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To  eliminate  much  of  this  difficulty,   equations  39  and  40  were 
developed.     A  similar  iteration  procedure  to  that  of  equations  47  was 
initially  used  to  investigate  the  new  equations,   and  it  was  found  that 
an  iteration  interval  of  0.  01  radians  gave  satisfactory  results  for    RQ 
between  0.  01  and  100. 

Equations  39  and  40  appeared  to  be  the  more  satisfactory  of 
the  two  sets  of  equations  for  programming  on  a  digital  computer,   since 
the  starting  solutions  of  numerical  methods  are  closely  related  to  the 
simple  iteration  procedure. 

A  Fortran  statement  for  the  solution  of  equations  39  and  40 
was  written  utilizing  a  subroutine  prepared  by  Mr.   K.   H.   Redner  of 
Computer  Sciences  Corporation.     The  subroutine  applies  the  method 
of  Plunge- Kutta- Gills  as  a  starting  solution,   and  the  Adams-Moulton 
predictor  -  corrector  formulas  as  a  continuing  solution.     The  facilities 
of  the  Computer  Center  at  the  California  Institute  of  Technology  were 
used  in  writing  the  program,   and  the  IBM  7090  computer     at  the  Jet 
Propulsion  Laboratory  was  used  for  the  actual  solution.     A  Fortran 
statement  is  included  as  Appendix  A. 

To  investigate  the  general  effect  of  changes  in  Y    and   n    , 
Stokes  drag  law  v/as  assumed  to  apply,   and  solutions  were  obtained 
for  six  values  of     F^o     between  0.  01  and  50  for  three  sets  of  assumed 
constants. 

M  Me    ■=■  JT     ,      Y=  I-*     ,      n-   2.o     , 

2)  Me     =nAF      ,      r  -    1-4-     ,       n  =  1.0       , 


3)         N\a  ■=.  J 2      .       r=  /-2      ,    n=  2.o 
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The  exhaust  conditions  of  a  typical  large  rocket  were  assumed 
as  initial  conditions  for  the  investigation  of  the  effect  of  allowing     -f 
to  vary  according  to  equation  23.     The  initial  conditions  chosen  were: 

M0  =    3.0      , 

7     Ibf  -  Sec 
/JLo    =   9.55XIO        ——r^ • 

ft       -   0.0.38     ^     , 
y     =    1.25"    , 

n     =  /.o     • 

Solutions  were  obtained  for    R0     of  0.  1  ,    1.0,   and  50  for  two  par- 
ticle sizes,  <J~    =  i  micron  and   <j~  =  2.  0  microns. 
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V.   RESULTS  AND  DISCUSSION 

The  results  of  this  investigation  are  plotted  on  Figs.    3  through 
10  for  0     to  70     of  flow  deflection  angle.     Flow  deflection  angle,    S    , 
was  chosen  as  a  more  meaningful  parameter  than    •&■    for  plotting  and 
analyzing  the  results  since  it  is  determined  solely  by  the  physical 
boundaries  of  the  flow.     The  fan  angle  and  deflection  angle  can  be  re- 
lated by  the  usual  Prandtl-Meyer  functions. 

Figure  2  is  a  graphical  representation  of  equation  6  giving     RQ 

as  a  function  of     To     and    <r    for  particles  having  a  density  of 

3 
210  lb/ft     (as  for  Al~Q-  ).     Initial  gas  characteristics  plotted  are  for 

air  at  100   F  and  1000   F  and  for  the  assumed  rocket  exhaust  conditions 
given  in  the  previous  section.      Ff0    is  not  greatly  affected  by  a  change 
in  initial  gas  conditions,   but  is  primarily  a  function  of   r^   and     <J~ 
For  particles  between  1.  0  u  and  4.  0  u  in  diameter  and  initially  dis- 
placed by  several  cm.   from  the  flow  boundary  (as  in  the  case  of  a 
rocket  exhaust),      R0   will  normally  be  greater  than  10. 

The  data  plotted  in  Figs.    3,   4,   and  5  is  based  on  the  assump- 
tions that   Mo  -  v2      ,    Y  =■  /.  4-    ,     n   -    2.  O   ,   and  Stokes  drag  law  ap- 
plies with  the  factor    -f  -  /.  O  .     por  these  conditions,  particles  with 
R.a  >  5  O      will  have  little  deviation  from  the  streamline,   and  slip 
velocities  will  be  low.     For     fR,0   <  O-  O]       there  will  be  large  tan- 
gential slip  and  very  little  deviation  from  its  original  flow  direction 
through  the  first  25     of  gas  deflection.     Radial  slip  is  small  compared 
to  tangential  slip,   being  a  maximum  of  .  164  for     f^0  —  o.  ol   at  about 
25     flow  deflection.     Near  this  point  the  maximum  relative  Mach 
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number  is  .  993  and  maximum  Reynolds  number  for  a  4.  0  ia  particle 
is  8.  38.     For     R0  -  j  o    the  maximum  radial  slip  is  less  than  one  per 
cent  while  the  maximum  Mach  number  is  only  0.  151  and  maximum 
Reynolds  number  for  4.  0  u  particle  is  2.  02.     The  Mach  numbers  and 
Reynolds  number  bath  drop  sharply  on  either  side  of  the  maximum 
values. 

The  assumed  viscosity  -  temperature  function,   equation  19 

£.-(!.)*,  (i9) 

/Wo  V  To   J 

is  a  well-established  approximation  that  applies  over  a  limited  tem- 
perature range  when  n   is  assigned  the  constant  value  appropriate  to 
this  range.     By  allowing    n    to  take  its  expected  limiting  values,   its 
effect  on  the  solution  can  be  observed.     Figures  6  and  7  show  curves 
for  which  Stokes  drag,      y  -  I.  4     ,   and   M0  ■=  v/E        v/here  assumed  for 
n  ■=.  I .  o  and    \r\  -  2. .  o       .     A  small  change  in  trajectory  path  and  tan- 
gential slip  velocity  and  a  negligible  change  in  radial  slip  velocity 
were  observed.     Thus,   little  error  is  introduced  by  the  selection  of 
an  average  value  of   n   over  an  extended  temperature  range. 

A  change  in    "2"    from  1.  4  to  1.2  does  have  a  fairly  large  ef- 
fect on  the  trajectory  path  which  is  due  principally  to  the  change  in 
the  gas  streamline.     The  tangential  slip  velocity  is  increased  also, 
which  tends  to  increase  the  path  deviation  from  the  streamline,   but 
this  is  a  small  effect.     Radial  slip  velocity  is  not  appreciably  affected. 

The  drag  correction  factor,    -f  ,   which  allows  a  departure  from 
Stokes  drag  law  to  compensate  for  slip  flow  or  free  molecular  flow 
effects,   was  proposed  by  Rannie  {ref.  3  )  as  an  analogue  to  a  similar 
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development  for  heat  transfer  coefficients.     Although  it  has  not  been 
verified  experimentally,   it  does  give  proper  results  for  the  extreme 


values  of     _Il   >>    /.  o         and    Ml    s /   \   o      »  and  should  give  a  good 

"e  ^<  e 

approximation  for  intermediate  values.     The  factor      — -        can  be 


shown  to  be  (ref.  7  ) 


Re  V  7T  r     Z<5~ 


where  $    is  the  mean  free  molecular  path  of  the  gas.     Thus,   if 

Mi    >>    I  •  O  ,   free  molecular  flow  exists,   and  if      Ml    <<  /.  o 

Re  Re 

continuum  flow  exists.     The  area  between  these  limits  is  the  slip 

flow  regime. 

From  equation  27  we  have,   for        Y=  '•  2    and     n  -  1.  o    , 


Ms  1     f       c°s   °< 

(48) 


R «,  Vcos(xe-K^) 


For  air  or  rocket  exhaust  gases  at  moderate  temperatures,   and  for 
particles  the  order  of  a  micron,     K,    will  lie  between  0.  03  and  0.  3  . 
Thus,   a  particle  of  this  general  9ize,  passing  through  a  Prandtl- 
Meyer  expansion,   will  start  in  the  slip  flow  regime  and  pass  into  free 
molecular  flow  at  large  expansion  angles  as  the  enclosed  cosine  func- 
tion in  equation  48  becomes  large  compared  to  unity. 

Figures  3,    9,    and  10  show  the  effect  of    -C    on  the  trajectory 
paths  and  slip  velocities  for  1.  0  micron  and  4.  0  micron  diameter 
particles.     The  deviation  from  the  uncorrected  Stokes  drag  law  tra- 
jectories is  large  for  particles  as  small  as  4.  0  micron  diameter, 
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and  increases  markedly  as  the  particle  3ize  is  reduced  to  1.  0  micron. 

In  the  case  of  a  1.  0  micron  diameter  particle,   an     R0  -  5  O   does  not 
closely  follow  the  streamline  as  it  did  for  Stokes  drag,   and      P?6  -  o.  ) 
now  shows  little  deflection  from  its  original  flow  direction.     For 
RQ  —  O.  \      a  maximum  relative  Mach  number  of  0.  924  and  maximum 
Reynolds  number  of  6.  73  is  reached  for  a  4.  0  micron  particle.     This 
compares  to  maximums  of  0.  853  and  6.  51  respectively  for  Stokes 
drag. 

The  maximum  values  of  Reynolds  numbers  obtained  are  greater 
than  the  limits  set  on  equation  23  for  the  determination  of   -f     .     A 
Reynolds  number  of  8.  0  would  change  the  equation  to 

-£-  =    0.65+    4.0  577    —    ' 
T  he  e 

(ref.  3  ),  which  would  make  a  23  per  cent  error  in  -f    for  a  Mach  num- 
ber of  0.  95  .     While  this  is  an  appreciable  error,   it  is  effective  over 
a  small  portion  of  the  total  expansion  and  diminishes  rapidly  on  both 
sides  of  maximum    Re    ,    so  it  should  not  greatly  affect  the  overall 
results. 
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VI-  CONCLUDING  REMARKS 

In  the  flow  regime  where  Stokes  drag  law  applies,   the  tra- 
jectories of  spherical  particles  passing  through  a  Prandtl- Meyer  ex- 
pansion fan  can  be  specified  by  a  single  dimensionless  parameter,  R    , 
which  is  a  function  of  initial  gas  conditions,  initial  radius,   and  particle 
diameter.     If  the  Reynolds  number  is  greater  than  about  0.  5,   the  tra- 
jectories are  also  a  function  of    Re  .       .hen  the  mean  free  path  of  the 
gas  molecules  is  of  the  same  order  or  large  compared  to  the  particle 
diameter,   Stokes  drag  law  must  be  modified  to  correct  for  the  discon- 
tinuous nature  of  the  flow,   and  the  trajectories  then  become  functions 
of  three  parameters,   R     ,   Re  ,   and    M  /Re  . 

It  was  found  that  the  temperature  dependence  of  viscosity  need 
not  be  precisely  known  to  ensure  satisfactory  results,   but  changes  in 
the  specific  heat  ratio,   y  »  c^°  appreciably  r-ffect  the  trajectory  paths. 

The  radial  velocity  lags  are  in  general  very  small  compa:. 
with  the  tangential  velocity  lags,   and  can,   in  many  cases,   be  neglected 
completely. 
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APPENDIX    A 
FORTRAN    STATEMENT 

CIKENSICN  Y(1C),  YCCTI3),  C(1C),  C(10) 

CIKENSICN  YUO),  YPCTI3),  C(1C),  C(1C) 
5t0  FCPfAl  (7F10.4) 
51C  FCRfAT  (1CH1  GAKKA  =  E13.6,  13H      MZERL  =  E13.6, 

19H      N  =  F5.2,  13h      RZERC  =  E13.6,  13H      OELTA  =  E13.6//) 
1C  READ  INPUT  TAPE  5,500,  CGAPPA,  CMZERO,  CN,  RZERO.DELT,  CMRE,  CO 

WRITE  CUTPUT  TAPE  6,  510,  CGAPKA,  CKZERO,  CN,  RZERO,  DELT 

Yd)  =  1.0 

Y(2)  =  1.0 

Y(3)  =  O.C 

T  =  0.0 

KAY  =  1 

CLAMCA  =  SCRTFI ( CGAfMA- 1 .0 ) / ( CCAKPA+ I . 0 ) ) 

E  =  (1.0/(CLAPCA«»2) l-CN+l.C 

ARC  =  SORTF  1 1.0+CLAMDA*CLAKDA«(CMZERO»CMZERC-1.0) ) 

ALPHA  =  ACCSFI l.C/ARG) 

C(l)    =     (  l.C/CLAPCA)-CLAPOA 

C(2)  =  C(  11/CLAKCA 

CO)  =  CCSF  (ALPHA) 

C(4)  =  CLAfCA»(C«3)«»(B-l.O) )»RZER0 

C(5)  =  2.0/CLAKDA 

C(fc)  =  CMRE  /  CD 

WRITE  CUTPUT  TAPE  6  ,  500  ,CL APCA , B , ARG, ALPHA 

WRITE  CUTPUT  TAPE  6  ,  500  ,C  I  1 )  ,  C  (  2  )  ,C  I  3  )  ,C  <  4  >  ,C  (  5  )  ,  C(6) 

CALL  CEC(K,3,T,Y,YDCT,.Cl,5.CE-7) 

GC  TO  t  1CO,20C,2CO,30C) ,K 
ICC  C(l)  =  CLAPCA.T+ALPHA 

C(2)  =  TANF  (C(  1  )  ) 

C(3)  =  CCSFICI  1  )  ) 

CIM  =  (  l.C/(r(3)»»B)) 

C(5)  =  l.C/(  1.C  +  CI6)  •  l(C(3)/C(3))  ••(B-1.0))) 

YCCTI  1)  =  Y( 3) 

YCCT13)  =  (C(5))  »Y[3)  *D(2)  -  C  (  2  >  •  Y  (  I  )  «0  (  2  )  *D  (  2  ) 
l-(C(4)»Y(3)/( (Y(1)»«3)*YI2)*Y(2) ) )»D(4)  •  0(5) 

YDCT12)  =  -C(  1)«Y(2)»C(2)+C(A)»0(A)«( 1.0/1 (Yl 1)»»3)»Y(2)) 
1-1.0/(Y(1)*Y( 1) ) )  *C(5) 

IF  (KAY-2)  150,  160,  16C 
160  CALL  HEG2 
15C  WRITE  CUTPUT  TAPE  6,  500,  C(l),  D(2),  D(3),  0(4),DI5) 

WRITE  CUTPUT  TAPE  6,  60C 
6CC  FCfiMAT  I//6H  THETA.9H 
11AH  Z,  16H 

214F  R/RC,  13H 

KAY  =  2 

GC  TC  160 
2CC  Y  (4  )  =  Y( 1 )»Y(2) 

Y(5)  =  Y(M«(  1.0-CLAMDA«YDCT(1  )/(  Y(  I  J»C<2)  )  ) 

Y(6)  =  (  (C(3)/C(3)  )«MCI2)  +  1.0)  )/Y(  I) 

Y(<5)  =  (  I  I.C-YI4)  )*«2)*l( ( 1.C-YI5) )«D(2)/CLAMDA)»»2) 

Y(7)  =  SCRTFI  Y(S))  ) 

Y(8)  =  Y(7)«( (C(3)/C(3) )•* (B-1.0  )  ) 

WRITE  CUTPUT  TAPE  6,  61C,  T,  ( Y I  JACK ) , JACK  =  1,8) 
610  FCRMAT (/F7.3,8E14.6) 
2C5  IFU-3. 0)210, 210,220 
210  CALL  CEQl 

GC  TC  205 
220  GC  TC  10 

CIPENSICN  Y(1C),  YDCT13),  C(IC),  D(IO) 

3CC  WRITE  CUTPUT  TAPE  6,  620 

620  FCRPAI  (22H  ERRCR  RETURN  FRCP  CEQ) 

GC  TC  10 

ENCI  1,0, 0,0, 0,0, CO, 0,1, 0,0, 0,0,0) 


1AH 

Y, 

V/VS, 1AH 

U/US, 

MS,  13H 

RE//) 
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